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The quantum kinetics of photons in hot QED plasmas is studied directly in real time 
^J , by implementing the dynamical renormalization group. In contrast to conventional ap- 

proach, the dynamical renormalization group method consistently includes off-shell (en- 
ergy non-conserving) effects and accounts for time-dependent collisional kernel. To lowest 
VQ ' order we find that in the relaxation time approximation the semihard photon distribution 

function relaxes with a power law. 



£S| ' The quantum kinetic description to study Abelian and non-Abelian plasmas in 

extreme environments is of fundamental importance in the understanding of the for- 
2^: ' mation and evolution of a novel phase of matter, the quark-gluon plasma (QGP), 

£\i ■ expected to be produced in ultrarelativistic heavy ion experiments. The typical ap- 

0*N proach to derive quantum kinetic equations begins by introducing a Wigner trans- 

form of a nonequilibrium Green's function at two different space-time pointsEJ and 
often requires a quasiparticle approximations The rationale behind the Wigner 
transform is the assumption of a wide separation between the microscopic (fast) and 
the kinetic (slow) time scales, typically justified in weakly coupled theories. Nev- 
' , ertheless, the quasiparticle picture, closely related to the assumption of completed 

qj ' collisions, is less warranted in ultrarelativistic heavy ion collisions, as estimates 

rl^ . based on energy deposited in ±he central collision region at BNL RHIC energies 

^/s ~ 200 A GeV suggest thatEl the lifetime of a deconfined phase of quark-gluon 
plasma is of order 10 — 20 fm/c with an overall freeze-out time of order 100 fm/c. 
We consider a hot, neutral QED plasma which is prepared at time t — to with 



fermions in thermal equilibrium at temperature T but photons out of equilibrium. 
The initial nonequilibrium photon distribution functions are given by n.k(*o)- The 
question that we want to ask is how does nu.(t) evolve in time? To make the 
calculation tractable we will focus on semihard photons of momentum eT <C k <C T 
(e <C 1 is the electromagnetic coupling constant), as the photon self-energy for 
semihard photon is dominated by the hard thermal loopQ (HTL) contribution which 
is still perturbative. 

To find out the evolution of n\^(t) directly in real time with initial condition 
n k(£o) specified, we set up an initial value problem for nk(t). This can be achieved 
as follows. First we construct a (Heisenberg) number operator Nk(t) that counts 
semihard photons of momentum k in terms of the (transverse) photon field operator 
and its conjugate momentum, and find dNk(t)/dt by using the Heisenberg equations 
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of motion. Then we express dn%c(t)/dt = (dNk(t)/dt) in terms of the nonequilibrium 
expectation value of the photon and fermion fields in nonequilibrium field theory 
and compute perturbatively in powers of e using the nonequilibrium Feynman rules 
and real-time propagators. To lowest order we findu 

j t n k (t) = [l + n k (( )]rf(()-n k (t )r>((), (1) 

rJ W . f"^*fo) H( " :*"'-''■" ■ (2) 



Here T^ (t) is the time- dependent rates and lZ k (u>) has a physical interpretation 
in terms of the off-shell (i.e., energy non- conserving) photon production (absorp- 
tion) processes in the plasma: bremsstrahlung e~(e + ) <-> e~(e + )7 and annihila- 
tion e~ e + «-► 7. Since the fermions are in thermal equilibrium, the Kubo-Martin- 
Schwinger (KMS) condition TZ^(lu) = e" /T TZ^lu) holds. In the HTL approxi- 
mation a detailed analysisO reveals that TZ^{uj) is completely determined by the 
off-shell (energy non-conserving) processes e~(e + ) <-> e~(e + )j. 

It should be emphasized that in the infinite time limit, the resonance factor in 
Eq. (0) can be approximated by a delta function, viz., 

sin[(w - k)(t - t )] t^oo 

w 5(w-fc), 

7T((jJ — fc) 

which is exactly the assumption of completed collisions invoked in time-dependent 
perturbation theory leading to energy conservation and Fermi's golden rule (i.e., a 
time-independent rate). In this limit the above off-shell (energy non-conserving) 
processes cannot take place due to kinematics, and the kinetics of photons should 
arise from higher order contributions. 

For any finite time, however, Eq. (|J) can be solved by direct integration, we find 

n k (t) = n k fo) + [1 + n k (to)] I dt'T<(t')- n k (t ) / dt'T>(t'). (3) 

The integrals that appear in the above expression, 

Jt J -00 (w - k) 2 

are dominated, in the long time limit, by the regions of u> for which the denominator 
is resonant (i.e, 10 ss k). Using the HTL approximation for TZ^(u>), we find thatEl 

dt'T<(t') fe(t -^ )>>1 ^El { l n [2k (t - to)} + 7 - 1} + 0[l/k(t - t )}, (4) 

to 0K 

where 7 = 0.577 ... is the Euler's constant. 

Do we find the time evolution of n k (t)? No, because at large times the (loga- 
rithmic) secular term in Eq. (|4j) that grows in time will invalidate the perturbative 
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solution Eq. (ph. The next step of our kinetic approach is tojesum these secular 
terms by implementing the dynamical renormalization groupJM which introduces 
a renormalization of the distribution function via 

M*o) = Zkfo, r)nk(r), Z k (t ,T) = 1 + e 2 z[}\t ,T) + e A z^ \t , r) + • • • . 

The renormalization coefficients z£ (to,T) are chosen to cancel the secular diver- 
gence order by order at an arbitrary time r (see Wang et alB for details). The 
perturbative solution can be evolved to large times provided that r is chosen close 
to t. A change in the time scale r is compensated by a change of the nk("r) in such 
a manner that n k (t) does not depend on the arbitrary scale r. This independence 
leads to the dmnamical renormalization group equation, which consistently to order 
e 2 is given byQ 

^n k (r) = [1 + n k (r)] F< (r) - n k (r) F> (r) + 0{e A ). (5) 

Choosing r to coincide with t in Eq. (H), we obtain the quantum kinetic equation 

jn*{t) = [1 + n k (i)] r<(t) - n k (*) r>(t). (6) 

In the relaxation time approximation, which describes the approach to equilib- 
rium of a slightly off-equilibrium mode with initial distribution n k (io) = ns{k) + 
<5n k (£o) [^s(fc) is the Bose-Einstein thermal distribution] while all the other modes 
are in equilibrium, solving Eq. (o) we find a power law relaxation at large times to 
be given byu 

Sn k (t)=6n k (t a )[k(t-t )}- e2T2 / 6k2 for k(t - to) > 1. (7) 

To conclude, our approach to derive kinetic equations in hot plasmas is different 
from the one often used in the literature which involves the Wigner transform and 
the assumption of completed collisions. In particular, it reveals clearly the dynamics 
of off-shell (energy non-conserving) effects arising from the finite system lifetime. 
This aspect is of phenomenological importance in the study of the direct photon 
production from a quark-gluon plasma of finite lifetime. u 
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